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We present a theory of collective spin excitations in diluted-magnetic-semiconductor quantum 
wells in which local magnetic moments are coupled via a quasi-two-dimensional gas of electrons or 
holes. In the case of a ferromagnetic state with partly spin-polarized electrons, we find that the 
Goldstone collective mode has anomalous k 4 dispersion and that for symmetric quantum wells odd 
parity modes do not disperse at all. We discuss the gap in the collective excitation spectrum which 
appears when spin-orbit interactions are included. 



PACS numbers: 75.50.Pp, 75.30.Ds, 73.43.-f 



c3 



i 

-a 
c 

o 
o 



(N 
> 

m 

(N 

m 
O 
-i— > 



I 

C 

o 
o 



X 



I. INTRODUCTION 

In the emerging field of spin- or magneto-electronics^^ 
the role of the spin degree-of-freedom in the properties of 
electronic systems is exploited in the design of new func- 
tional devices. The recognition of this additional degree 
of freedom suggests possibilities for electrical manipula- 
tion beyond the tool-set of conventional electronics which 
is based entirely on coupling to the electronic charge. The 
effort to generate and manipulate spin-polarized carri- 
ers in a controllable environment, preferably in semicon- 
ductors, has triggered the discovery of carrier-induced 
ferromagnetismi^ in diluted magnetic semiconductors 
(DMSs), 5 In these systems a few percent of the cations 
in III-V or II- VI semiconductor compounds are randomly 
substituted by magnetic ions, usually Mn, which have lo- 
cal magnetic moments. The effective coupling between 
these local moments is mediated by free carriers in the 
host semiconductor compound (holes for p-doped ma- 
terials and electrons for n-doped ones) and can lead to 
ferromagnetic long-range order. Curie temperatures T c 
in excess of 100 K have been found in bulk (Ga,Mn)As 
systems/ " 

One approach to understand the magnetic and opti- 
cal properties of DMSs is based on a phenomenological 
model of the relevant low-energy degrees of freedomi&S 
In this picture, local S = 5/2 spins^ from Mn 2+ ions 
are exchange coupled to itinerant carriers of a metallic 
nature. In typical samples, the density of free carriers 
is much smaller than the Mn ion concentration. For n- 
doped materials, the exchange is due to ferromagnetic 
s-d coupling, while for p-doped ones it is due to antifcr- 
romagnetic p-d coupling, as illustrated schematically in 
Fig. ^ hi both cases, the free carriers are believed to 
mediate an effective ferromagnetic coupling between the 
Mn spins, which is typically stronger than the shorter- 
range antiferromagnetic direct exchange coupling present 
in undoped systems. 

The reliability of this phenomenological approach has 
been tested by comparing theoretical predictions with 
experimental findings. The tendency towards ferromag- 
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FIG. 1: Schematic representation of the exchange coupling 
between itinerant-carrier and localized magnetic-impurity 
spins in n-doped and p-doped DMSs. When the local mo- 
ments are parallel to each other and the band system is spin- 
polarized, the exchange energy can be minimized for either 
ferromagnetic or antiferromagnetic interactions. 



netic order and trends in the observed T c 's, domain struc- 
ture properties, the anomalous Hall effect, and magneto- 
optical properties, have been successfully described by 
treating this phenomenological model in a mean-field 
approximation (MFT), 13 i 14 i 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 which is 
analogous to the Weiss mean-field approach for lattice 
spin models. In the mean-field theory the local Mn ions 
are treated as independent but subject to an effective 
magnetic field which originates from their exchange in- 
teractions with spin-polarized free carriers. Similarly, 
the itinerant-carrier system sees an effective field propor- 
tional to the Mn density and polarization. This picture 
does not account, however, for correlations between Mn 
spin configurations and the itinerant carrier state which 
reduce the energy cost of local-moment spin fluctuations 
that have slow spatial variations. As a consequence, 
MFT systematically overestimates the Curie tempera- 
ture, a problem which is severe for systems with re- 
duced dimensionality^ including the quantum well sys- 
tems that will be discussed here. 

One prediction that follows from the phenomenologi- 
cal model is that the system's collective excitations in- 



2 



volve correlated dynamics of local moment and itiner- 
ant spins. In the case of bulk DMS systems, we have 
predicted two branches of collective spin waves and dis- 
cussed their properties^ as well as their impact on lim- 
iting the Curie temperature,— This analysis of collec- 
tive excitations requires a theoretical description beyond 
MFT, which neglects correlations, and beyond the famil- 
iar Ruderman-Kittel-Kasuya-Yoshida (RKKY) theory of 
pair-wise carrier-mediated interactions, which fails for 
the systems under consideration because it assumes a 
carrier-band spin splitting that is small compared to the 
Fermi energy. This assumption is not typically satisfied 
in doped DMS systems^ in part because the itinerant- 
carrier concentration is usually much smaller than the 
Mn impurity density^ Moreover, the RKKY picture also 
assumes an instantaneous static interaction between the 
magnetic Mn ions, neglecting the retarded character of 
the itinerant-carrier response that mediates the interac- 
tions. 

An indication 2 ^ that the spin excitations of doped 
DMS systems have collective local-moment and car- 
rier character, even in paramagnetic systems, has been 
provided by recent electron paramagnetic resonance 
experiments 29 in n-doped DMS quantum wells. The aim 
of the present paper is to extend the previous theoreti- 
cal work to describe the full dispersion of all collective 
spin excitations in quantum wells, their dependence on 
the magnetic-ion doping concentration and profile, and 
on the free-carrier density. This is a first step toward the 
theoretical study of quantum and thermal fluctuations in 
the magnetism of nanostructured DMSs which are start- 
ing to receive increased attention, partially because of 
the possibility of quantum confinement control of mag- 
netic properties as in the recent experimental study of 
a DMS quantum well in Ref. IsoL Quantum confinement 
is expected to drastically affect the magnetic properties 
of nanostructures, 31 ! 32 ! 33 ! 34 ! 35 ! 36 ! 37 ! 38 ! 39 ! 40 In Ref. |H only 
the long-wavelength limit of the lowest spin-wave branch, 
the mode that electron paramagnetic resonance probes, 
was considered. 

The article is organized as follows. In Sec. [H] we 
develop the theoretical tools necessary to address col- 
lective excitations in doped DMS quantum wells in a 
general way. After introducing the many-body quan- 
tum Hamiltonian (Sec. Ill A(l we derive an effective ac- 
tion (Sec. Ill Bll that leads to an independent spin-wave 
theory for low temperatures in multi-subband quantum 
wells fSec. III(J|l . Sec. IIIII is dedicated to the evaluation 
and discussion of collective spin excitations, concentrat- 
ing on the case in which a single electronic subband is 
occupied and subband mixing is negligible. We find that 
odd-parity collective modes of doped quantum- well DMS 
systems are dispersionless in this limit, and that the low- 
est energy Goldstone collective mode of ferromagnetic 
systems has anomalous k 4 dispersion when the quantum- 
well carrier system is not half-metallic (i.e. when the car- 
riers are not fully spin-polarized). Results for dilute and 
moderate Mn doping are shown in Sees. IIIlO and lTil Dl 



respectively. The role of spin-orbit coupling, which gives 
rise to magnetic anisotropy and creates a gap in the ex- 
citation spectrum of a ferromagnetic system, is discussed 
in Sec. IIII El A summary and discussion of our results is 
presented in Sec. IIVI 



II. DERIVATION OF THE THEORY 

A. Hamiltonian 

We consider a symmetric quantum well of uniform 
width d that confines the motion of itinerant carriers 
in the z-direction (see Fig. [5J we later comment the 
case of asymmetric quantum wells). The carriers move 
freely in the x-y plane, occupying one or several trans- 
verse modes or subbands. The quantum- well geometry 
makes it convenient to split the three-dimensional (3D) 
spatial coordinate into (r, z), where r corresponds to the 
two-dimensional (2D) x-y-projection. The field opera- 
tor for the itinerant carriers ^(r, z) can be written as 
= El'=i^(r)Xm(4 where m = 1,...,M la- 
bels the subband number, Xm(z) = y/2/d sin(m7rz/<i) is 
the real wave function for subband m, which satisfy the 
orthonormality condition f^dz Xm{z)Xm'{z) = 8 m , m >, 

and Vvn( r ) is a spinor with components Vv,m(r)- As in- 
dicated above we will adopt particle-in-a-box wavefunc- 
tions for explicit calculations, although this approxima- 
tion plays no critical role in our theory and can easily 
be relaxed. The transverse wave function Xm{z) degree- 
of-freedom will later be taken to be frozen in its ground 
state; this is normally a good approximation except in 
wide quantum wells. The in-plane degrees of freedom 
are described by fluctuating spinor fields 'ip a ;m( r )- The 
magnetic impurities are randomly distributed within the 
quantum well. 

The fact that the Mn density in typical quantum well 
systems is very much larger than the carrier density sug- 
gests the replacement of the random distribution of local 
Mn magnetic moments by a continuous density Nmu(z), 
thereby neglecting disorder in the Mn ion locations! 41 
This leaves us with a situation in which a growth di- 
rection degree of freedom exists for the local moment 
spins, but not for the quantum well electrons. It is the 
quasi-3D character of local moments that are coupled 
together by quasi-2D electrons that is responsible for un- 
usual aspects of the collective excitation spectrum that 
we will discuss later. The Debye-like continuum approx- 
imation we use for the Mn ion density distribution will, 
of course, fail for modes that involve either in-plane or 
growth-direction spatial variation on a scale shorter than 
the distance between Mn ions, as we discuss later. The 
dependence of Nyi n (z) on the growth-direction coordi- 
nate z allows for the possibility of a non-uniform doping 
profile in the quantum well. The two-component spinors 
we use for the quantum-well electron fields restrict our 
attention to circumstances in which the electric subbands 
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occur in pairs with identical orbital wavefunctions, i.e. to 
hole quantum-well subbands with small heavy-light mix- 
ing or to electron subbands. Generalizing to arbitrarily 
spin-orbit coupled systems considerably complicates the 
notation we use below and, in the case of an external 
field, complicates the theory considerably because of the 
interplay between orbital and Zeeman coupling. We will 
for the most part restrict our attention to n-doped quan- 
tum wells. 
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FIG. 2: Sketch of a DMS quantum well. Itinerant carriers 
move freely in the x-y plane, occupying subbands \n ( z ) due to 
quantum confinement along the z-axis (intersubband energy- 
gap Ad). The magnetic-ion doping profile is represented by a 
continuous Mn-density distribution Nmu(z). 

The total Hamiltonian H consist of four terms: H = 
Hkin + Hz + -Hex + Hy>. In the presence of a magnetic 
field B = V x A, the kinetic term for carriers of charge 
e reads 



Aid, 



d r / dz x 



*t(r, z) 



2m 



-V 2 + V(z) - fx 



tf(r, z) 



r M 



2m" 



(1) 



Vv,m(r), 



where V( r j = V( r ) — (ie/hc)A, V(z) is the quantum-well 
confining potential, and fi' m = /i — e m is the effective 
subband-dependent chemical potential of the quasi-2D 
carrier gas, with e TO the subband quantization energy. 
Here, we assumed a parabolic dispersion for the free carri- 
ers, with an effective mass m* . This is well justified for n- 
doped systems, which have s-band conduction electrons. 
For hole doped systems, which have p-band valence car- 
riers, this approximation is often useful for qualitative 
discussions. The Zeeman term is 



H z = A*bB 



dz [g e s(r,z) +g Mn S(r, z)] 



(2) 



where /iB > is the Bohr magneton, 
s(r,z) = -¥{r,z) t $(r,z) 



(3) 



E Yx m {z)x m '{z) ^ m (r)^^Vv, m ,(r) 



m,m' cr.cr' 



is the quantum-well carrier spin density (with Pauli ma- 
trix vector r), and S(r, z) is the spin density of the Mn 
subsystem. The coupling between the carrier spins and 
the local Mn spins is described by 



(4) 



#cx = J cx J d 2 r J dz S(r, z) -s(r,z), 

where J GX < corresponds to ferromagnetic and J cx > 
to antiferromagnetic coupling (i.e. to n- and p-doped 
host semiconductors, respectively). In symmetric quan- 
tum wells spin-orbit interactions are described by the 
Drcsselhaus Hamiltonian^ 

H D = 7 J d 2 r^ dz¥{r,z) k 2 z {-T x k x + T y k y )^{Y,z) 



M 
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(5) 



where (k 2 z )„ hm > = J Q dz Xm{z) k 2 z Xm'(z) = {k 2 n )8 m , m , 
with (fc^) = (mir/d) 2 when particle-in-a-box orbitals are 
used. The above spin-orbit Hamiltonian H^ leads to a z- 
oriented magnetic easy-axis as we have shown in earlier 
work2& and discuss later. For comments on spin-orbit 
coupling in the case of asymmetric quantum wells see 

sec. mrfei 



B. Effective action 

In analogy to our earlier workS4 on bulk DMS ferro- 
magnets, we want to describe elementary spin excita- 
tions in the DMS quantum well in a language where the 
itinerant-carrier degrees of freedom are integrated out. 
This leads to a retarded free-carrier mediated interac- 
tion between the Mn-ion S = 5/2 spins. We are in- 
terested in small spin fluctuations about the mean-field 
magnetic state. The ground state of experimental doped 
quantum well DMS systems have sometimes been found 
to be ferromagnetic; 30 ^ 1 ! 32 ^ 3 ! 34 ^ 5 ! 36 ^ 7 ! 38 ^ 9 ! 40 and some- 
times exhibit complex spin-glass behavior. It is quite pos- 
sible that the complex spin-glass states that sometimes 
occur are due to disorder effects that are not essential 
and can in principle be avoided, due for example to in- 
homogeneities in the Mn ion distribution, substitutional 
Mn ions, or other defects. In any event, the theory we 
discuss assumes a mean-field state in which all Mn ions 
are aligned. When this simple state is not the ground 
state of the system, or when we want to describe the col- 
lective excitations of a system that is above its ferromag- 
netic transition temperature, our theory will apply only 
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in an external magnetic field that is strong enough to 
achieve substantial Mn ion spin polarization. We choose 
this field to be oriented in the z direction: B = (0, 0, B). 
Because <?Mn > 0, the Mn spins then tend to align along 
the opposite direction. In the case of ferromagnets with 
anisotropy, we choose the z direction to be along an easy- 
axis. 

It is convenient to represent the S — 5/2 spins 
by Holstein-Primakoff (HP) 44 bosons. For small fluc- 
tuations around the mean-field state the spin density 
S(r, z) is approximated by S + « Coy / 2NMn(z)S, S~ w 
ujyj2Nyin{z)S, and S z — ujlo — Nyi n {z)S, where the com- 
plex variables Co, lo are boson creation and annihilation 



operators that become bosonic coherent state labels in 
the path-integral formalism we employ. The vacuum with 
no HP bosons corresponds to full (negative) polarization 
of the Mn system, while the creation of a HP boson de- 
scribes an increase in the total Mn spin by one unit. The 
partition function Z of the compound system is calcu- 
lated using a coherent-state path integral representation 

Z = J V(^)V(Qu) e~ dT ( 6 ) 

with T>(ipil>) = Ufyiipi) . . . T>(ipMipM) and the La- 
grangian 



L(tpip,Lou) = d 2 r i> m (T,T)d T il) m {r,T) + / dz w(r, z, r)d T Lo(r, z, r) 

J lm=l J " 



(7) 



r 



where the Grassmann numbers ip — (ipi, . . . , 4>m ) an d 
if) = ("01, . . . , iPm) describe fermion (itinerant carrier) 
fluctuations within each of the M subbands. 

Since the Hamiltonian is bilinear in the fermionic fields, 



we can integrate them out and arrive at a represen- 
tation for the bosonic partition function of the form 
Z = JT)(u>u))exp(—S e s[u>u>]), with an effective action 



S e s_[ujuj] = I dr / d 2 r / dz [u)(r, z, r)d T Lo(r, z, r) + gMnpsB • S(ww)] — ln[det G~ (Qui)]. 
Jo J Jo 



(8) 



r 



The total kernel G _1 ((I'ijj) can be split into a mean-field 
part, which does not depend on the bosonic (Mn spin ex- 
citation) fields Co and u>, and a fluctuating part by writing 
G- 1 {luuj) = GI 



-l 

'MF 



(^MFj 



SG (uiu>), with 

h 2 



3r- 



2m* 



(9) 



1 



^ (SeA*B-B<W, 

+ j( k m) {-T X kx + Tyky) 5; 



>)t z 



y y ) u m,m' j 



6G- 



Je: 



(UJLO 



(z) X 



(10) 



\J2Nmti{z)S (CO T +10 T + ) + LOLOT Z 



The exchange coupling contributes to the conduction 
band spin-splitting in Gmf through the mean-field in- 
teraction A m>m ' = J CK {NMn)m,rn'S, where (7V Mn )m,m' = 

f Q dz Xm(z)Xm'(z)NMn(z) 45 We recognize here that 
coupling to a Mn-spin system with an inhomogeneous 
doping profile Nyi n {z) leads to mixing of the quantum- 
well subbands. These intersubband interactions are 
present at the mean-field level, as seen in Eq. 10, but 



also appear in the term which expresses the coupling be- 
tween carriers and local moment fluctuations, Eq. (|10fl . 



Quantum-well subband decoupling and 
independent spin-wave theory 



The above picture simplifies considerably when 
quantum- well subband mixing is negligible^ i.e. when 
the energy gap A^ (see Fig. [5} is much larger than 
the spin-splitting energies A TO)m >. This regime can be 
reached either by narrowing the quantum well (A^ — > 
oo), or by diluting the Mn doping (A m>m / — > 0), indepen- 
dently of the number of occupied subbands M , which is 
controlled by the carrier density. In this limit we arrive 
at a Green's function that is diagonal in subband space 
with G- 1 = (GzL) m + 5G 



-l 

* m ' 



MFi™ 



Or 



2m* 



(11) 



+ -zideimB - A m )T z 

+ 7(fcm) (-T X k X +Tyky) , 



5 



5G-\Qu) = J ~f dz X 2 m (z)x (12) 

y/2NMn(z)S (U> T~ + LO T + ) + U}UJT Z . 

Here A m = A m m is the subband-dependent exchange 
contribution to the itinerant carrier mean-field splitting. 
This approximation leads to a convenient subband sepa- 
ration lnfdetG- 1 ] = Em^Mdet^m 1 ] in E q- ©■ 

Expanding ln^etG" 1 ] up to second order in 5G" 1 , 
lntdetG" 1 ] = tr[ln(GMF) m ] + tr[(G MF ) m '5G m 1 ] - 
(l/2)tr[(G M F)m'5G m 1 (GMF)m5G^ :L ] + . . ., and collecting 
all contributions up to quadratic order in ui and u) we ar- 
rive at an independent spin-wave theory where the spin 
excitations are treated as non- interacting HP bosons. 
This is a good approximation for temperatures well be- 



low the maximum of the ferromagnetic transition tem- 
perature and/or the temperature defined by Zeeman cou- 
pling to the external field, in which case spin excitation 
amplitudes are small. Fourier transforming the resulting 
spin-wave action (keeping z in real space and defining 
bosonic Matsubara frequencies v n ) we obtain 

5cffM = I J 'W? Jo"" H dZ ' X 

w(k, z, j/„)Z?" 1 (k, k', z, z', i/„)w(k', z', u n ).(13) 

The kernel of the quadratic action (|13|) is the inverse of 
the spin-wave propagator Z?(k, k', z, z', u n ) and is given 
by 



A I 



-Wn + ffMnMS- 8 - ~7T Yl (™™ ~ ™m)Xm( z ) 



J 2 

CX ■ 

2 



m— 1 

M 



5(z - z')5{k - k') 



J --5 V^M„(z)iVM„(.') £ X 2 m {z)x 2 m {z')Y. fieLa \ *% a '*(k) ^'(k')- (14) 
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We arrive at Eq. I|14|l by introducing a wave- function rep- 
resentation of the mean- field Green's function of Eq. i|ll|) 

(G MF ) m (r ,T,V n ) = }^ _ , (15) 

where </>^ a (r) are the 2D mean- field itinerant carrier 
eigenstates for spin a and subband m, with energy Q = 

e a + (a/2)(g e /j, B B - A m )^l + ie s n ° a /(g e fi B B - A m ) - 

/j,' m mL The index a accounts for quantum- numbers asso- 
ciated with orbital motion, with kinetic and spin-orbit 
energies e a and e^ Q = \(l,m,a\'yh%(-T x k x + T y k y )\ f 
, m, a) \ 2 / (g e ^B — A m ), respectively. To first order in 
spin-orbit coupling (e^° a -C (g e ^BB — A m )), the eigenen- 
ergy is e£, a = e Q + (a'/2)(g e fj, B B ~ A m ) + ae™ a - fj,' m . 
In Eq. {T2J /(e) is the Fermi distribution, is the 2D 
mean-field itinerant carrier spin density for the m,a sub- 
band obtained by summing over occupied states, and we 
have introduced Fourier-transform factors defined by 

$r'(k) = J d 2 r exp(zk ■ r) 4,« (16) 

The first line on the r.h.s. of Eq. I|14l) is local in 
space. It represents the mean-field expression for the 
exchange field that the Mn spins experience. The sec- 
ond line is nonlocal in space and describes correlation 
effects that occur because of the (space- and m-subband- 
dependent) response of the quantum-well carriers to Mn 
spin orientations. We note that D~ x is not a function of 



z — z' only, because of the absence of translational sym- 
metry in the z-direction. Additionally, a 2D Debye cutoff 
fen = 47riVMn /{N/d) ensures that our continuum approx- 
imation has the correct number of magnetic-impurity de- 
grees of freedom. [Here, ./VMn = Jq dzNMn(z) / d — and 
N is the number of growth-direction modes included in 
the theory as we explain below. We associate N with the 
mean number of Mn ions encountered on crossing the 

1 /3 

quantum well; e.g. for isotropic doping N ~ JV Mn d.] 

We comment now on the factors (k) which are 
trivial for the plane- wave functions of field- free systems. 
They are included to allow us to simply account for the 
consequences of orbital coupling of itinerant carriers to 
magnetic fields. The index a includes both the Landau 
levels index and the gauge-dependent index for states 
within a Landau level. At zero field the mean-field eigen- 
states </>J^ a (r) are plane waves with momentum a = q, 
kinetic energy e q = h 2 q 2 /(2m*) and spin-orbit energy 
4°q = l 2 (kl l ) 2 q 2 /(gefiBB~A m ). In either case Eq. 
is diagonal in the in-plane momentum k and we denote 
its diagonal elements by L> _1 (k, z, z', iv n ). The case of 
a uniform magnetic field was considered in Ref. |2^ and 
will not be discussed further in this paper. 



III. ELEMENTARY SPIN EXCITATIONS 

Even the bulk-like epitaxially grown thin film samples 
studied in typical experiments do not contain a very large 
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number of occupied 2D carrier subbands. Our formal- 
ism could in principle be used to calculate the collective 
modes of thin films, taking account of the variation in Mn 
density across the film, although the approach becomes 
numerically cumbersome when more than a few sub- 
bands are occupied. It is likely, though, that greater in- 
sight into vertical inhomogeneity effects in thin-films can 
be obtained with more approximate approaches^ We 
limit the discussion here to true quantum-well samples in 
which a single subband is occupied (M = 1) and subband 
mixing can be neglected. For definiteness we concentrate 
on the case of constant Mn density NMn(z) = Nmu- Sim- 
ilarly we assume that the state we are studying is an or- 
dered ferromagnetic phase with the z-direction magnetic 
easy-axis that is favored by (weak) spin-orbit coupling. 
The situation in which the magnetization direction has 
been reoriented by an external magnetic field is readily 
included in the formalism as explained above. General- 
izations to the cases of multiple subbands, and inhomo- 
geneous Mn density are straightforward. 

Collective spin excitation dispersion f2(k) branches are 
located by finding the frequencies at which the determi- 



nant, det[£> 1 (k, 



/„ = f2)], of the quadratic ac- 



tion kernel in Eq. (|13|) vanishes. The continuum of 
spin-flip particle-hole (Stoner) excitations is located by 
identifying the k-dependent frequency range over which 
Im[.D _1 (k, z, z' , w-n)] is non-zero after the analytic con- 
tinuation iv n — ► O + i0 + . 



Stoner continuum 



We start by evaluating the continuum of Stoner exci- 
tations introduced above. They correspond to flipping 
a single spin in the itinerant carrier subsystem and typ- 
ically have relatively large energies of the order of the 
itinerant carrier mean-field spin splitting A = J cx NmuS. 
The continuum is obtained by determining the condi- 
tions for Im[D^ 1 (k,z 1 z',n + i0+)] ^ 0. With this aim 
it is convenient to define the dimensionless carrier spin- 



polarization p = (n 



T — ri.l 



)/(r 



) and the Fermi energy 



of the majority-spin carrier band e s = pi + |A|/2, where 
pi is the effective chemical potential of the 2D carrier gas. 
For half-metallic carriers (\p\ = 1, e s < |A|, see Appendix 
Ia"|) these excitations carry spin S z = ±1, depending on 
the sign of A (i.e. on whether the coupling between 
carriers and Mn ions is ferromagnetic or antiferromag- 
netic). On the other hand, for partly polarized carriers 
(\p\ < 1, e s > |A|, see Appendix Ia")| excitations carrying 
both S z = 1 and S z = — 1 contribute, independent of the 
sign of A. In the absence of spin-orbit coupling (7 = 0), 
one finds a continuum of excitations with dispersion lying 
between the curves — A — sign[A]ei c ±2 v / e^ek for e s < |A|, 
and also between — A + sign[A]e k ± 2yJ (e s — |A|)e k for 
e s > |A|. For small spin-orbit coupling, the energy width 



approaching the width 
{1 2 {k 2 z ) 2 m*e s 



An 



h 2 \A\ 



(17) 



In the case of multi-subband quantum wells (M > 1) 
multiple continua arise, each of them associated with the 
corresponding subband by means of A m , (fc^J and p' m . 



B. Spin- wave modes 

The number of collective modes that appear in our the- 
ory depends on the doping concentration and the width of 
the quantum well. It is natural to choose the mean num- 
ber N of Mn ions along the z-direction as a dimensional 
cutoff for the representation of the inverse propagator 
.D -1 (k, z, z', f2). This motivates the choice of an appro- 
priate basis of N orthonormal excitation profiles {uj n (z)}, 
with < n < N — 1 and J dz LU n (z)uj n i(z) = 8 nn >, for 
expanding Eq. (|14f> . i.e., 

d d 

D-*,(k,n) = J dz J dz'u n (z)D- 1 (k,z,z',Q)u n ,(z')(18) 


We later solve det[£)~^,(k, f2)]jvxjv = for fl and 
obtain a set of TV + 1 solutions {fim(k)} with the 



mode profiles u>m(z) 
cients : 



[^,(k,0 ( n)] 



JO r (0 



(1))\NxNC 



(1) 



^„cl Z) w„(z). The coeffi- 

' c tv-i) are obtained from 
0. We combine this procedure 
with a Debye cut-off of the 2D wavevectors to get the 
correct number of magnetic degrees of freedom. This ap- 
proximate procedure, a silent partner of the continuum 
Mn density approximation that we use to avoid dealing 
with disorder, obviously breaks down to some degree for 
the shortest wavelength modes which must be sensitive 
to the discreteness of the magnetic degrees of freedom. 
The procedure should be accurate for longer wavelength 
modes, however, and we believe that it gives a good qual- 
itative description of the overall spectrum. We employ it 
without further comment in the rest of the paper. 

We assume that the magnetization direction of the Mn 
spins located at the borders of the quantum well is not 
fixed by an anisotropy field or magnetic coupling to an 
adjacent layer. Then we can use free-end boundary con- 
ditions for the magnetic excitations, duj n (z) / dz — at 
z = 0,d. This determines the choice of the basis func- 
tions 



U} n (z) = a n cos \ ) , O-r, 




for n = 
for n > 1 



(19) 



We now calculate the matrix elements -D „, (k, f2) us- 
ing Eqs. ifH Il -lfTQ" ! for M = 1 and constant JVm n in the 
absence of an external magnetic field. The quantum- 
well subband is defined by the wave function x( z ) = 
y/2/d sm(-rrz/d). We find 



of the Stoner continuum does not vanish at k = 0, instead D nn , (k, 51) 



|A|[Mi + Z(k,n)M a ] 



(20) 
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with the TV x iV-matrices 



Mi = 



/ 1 





-i/V§ 











1/2 





-1/2 





-l/x/2 





l 





-1/2 





-1/2 





1 











-1/2 





1 



M 2 = 



1 -1/V2 



-1/V2 1/2 





V 



(21) 



NxN 



(22) 



NxN 



free- 
,T _ 



In Eq. (|20fl we have defined the ratio of the 
carrier spin-density to the Mn spin density x s = \n< 
n^\/ (2NM n Sd), which typically satisfies x s <C 1, and the 
dimcnsionlcss integral 



i(k,n) 



q 



/(e T 



q+k 



(23) 



The last can be evaluated analytically for T — in the 
absence of spin-orbit coupling, as described briefly in Ap- 
pendices El and IbI 

The term containing Mi corresponds to the term pro- 
portional to J cx in Eq. (|14f> and describes the mean- 
field exchange interaction between Mn spins with free- 
carrier spins. The appearance of off-diagonal matrix el- 
ements, indicating a mixing of basis functions for the 
Mn spin excitations, is of geometric origin, determined 
by the projection Jq dz x 2 (z)w„(2)w n /(z). The nonlo- 
cal correlations are accounted for by the term contain- 
ing M2, which corresponds to the term proportional to 
J 2 X in Eq. (|14f> . Mixing appears here also, determined 

this time by [j Q dz X* {z)v n (z)][J dz x 2 ( z )^n'(z)}- 
The structure of the matrices in Eqs. I|21ll and 
(I22|l shows that basis functions with different par- 
ity {uj n (z)} do not mix. 50 This allows us to write 
the expanded kernel (|20|l as the (external) product 
of two matrices corresponding to even (+) and odd 
(-) modes: [D^,(k,n)] NxN = [D~ x n , (k, Q)}+ + xN+ ® 

[£)-^(k,0)]-_ xiv _, with N+ + N- = N. Spin 

modes obtained as solutions of det[D~^,(k,fl)]ffxN — 
can now be classified according to their parity 
by solving separately det[Z)~*,(k, Q)}^ +yN+ = and 

det[D~^, (k, ri)]~_ xjv _ = 0, respectively. This leads to 
N + + 1 even modes and N~ odd ones. Moreover, it is 
possible to see from Eq. (|2*U|l that [D~*, (k, _ N _ 
is independent of /(k, 51). This means that correlations 
between local moment and band configurations do not 
influence spin modes of odd parity. As a consequence 
these modes are dispersionlessr 0, as we see explicitly be- 
low. More interesting are the even modes for which cor- 
relation effects due to the coupling to itinerant carriers 
show up. 



In the following we apply the above formulation to the 
cases of dilute and moderate Mn doping in the limit 
of vanishing spin-orbit coupling (7 — > 0). We then 
(Sec. " 



IIII Ell comment on how these results are altered 
by a finite 7. 



C. Dilute Mn doping 

For illustration we start discussing the limiting case of 
dilute Mn doping or, equivalently, narrow quantum wells. 
This corresponds to very few Mn ions across the quan- 
tum well leading in our approach to a low-dimensional 
kernel [D~^,(k,£l)] NxN (Eq. J20J) with N of order one. 
In this situation we can easily approach the problem an- 
alytically. We choose for simplicity N + — N~ — 1 (i.e. 
N = 2). In this case the dispersion relations of even 
and odd spin modes (corresponding to ujq(z) and u>i(z), 
respectively; see Eq. Hiyf) ) are obtained by solving 



det[I^,(k,fi)]- 
det[£>-^(k,n)]- 



[l + 7(k,n)]a: B |A| 
af s |A|/2-n = 0. 



n = 0, 



(24) 
(25) 



Eq. I|25(l leads to a single odd mode with flat dispersion 
fl~ = 3c 8 |A|/2. For the calculation of the even modes 
we limit ourselves for now to the case of ferromagnetic 
coupling (A < 0) at T — 0; the difference between ferro- 
magnetic and antiferromagnetic cases is commented on 
later in Sec. IIVI We solve Eq. (|2~4"|l for long and short 
wavelengths (i.e. long and short range correlations), us- 
ing the expansions (|B2p and l|B3J) up to first order in et 
and 1/ek, respectively. This leads to two solutions: one 
A|, and one hard mode f2^ iff ~ |A|, 
For half-metallic carriers 
(\p\ = 1, e s < |A|, see Appendix 0) we obtain for the 
small-and large-momentum limit 



soft mode f2+ ft < x. 



where typically f2+ ft <§; fi^ iff . 



^ s + of t(k) 



zJAI 1 



o(V4). 



n+K(k) = (i + x s )\a\ 



1 



1 + Xs 



IAI 



(26) 
(27) 

^ + 0(el), 
(28) 



respectively. Correspondingly, for partly polarized carri- 
ers (\p\ < 1, e s > |A|, using the results summarized in 
Appendix lA"|l we find 



ft+ oft (k) = + O(4) 



xJA\ 1 



2e s 

"1ai; < k 

1 /2e< 



(29) 

O(l/4),(30) 



n+ lff (k) = (1 + 2 ; S )|A| + - ^g-lj £k + (4),(31) 

for small and large momenta, respectively. 

The branch flt^. corresponds to a gapless Goldstone- 
mode reflecting the spontaneous breaking of rotational 
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symmetry for the Mn spins subsystem, as expected 
generically in ferromagnets and found in bulk magnetic 
semiconductors^ At long wavelengths (ek — > 0) the dis- 
persion in bulk isotropic ferromagnets is proportional to 
the spin stiffness p divided by the magnetization M (i.e. 
£1 = (p/Ai)k 2 ). Similarly, in the adiabatic limit e s <C 
|A|, our long wavelength result for quantum wells lj26|l 
shows a spin stiffness due only to the increase in kinetic 
energy of the fully spin-polarized band when the spin 
orientation is spatial dependent, p = h 2 n^ / (4m*). The 
magnetization, with parallel contributions from Mn ions 
and itinerant carriers coupled ferromagnetically, reads 
M = SN Mn d + n l /2 = SN Mn d(l + x s ). However, un- 
like bulk systems the spin stiffness vanishes as e s — ► |A|, 
Eq. pfijl. and stays equal to zero for e s > |A|, Eq. 
This unusual feature should lead to some non-standard 
phenomenology in these ferromagnets, for example in 
the physics that controls domain wall widths and finite- 
temperature magnetization suppression. For short wave- 
lengths (ek — > oo), Eqs. l(77l) and i|5U)l. the excitation en- 
ergy fi,i ft tends to a mean- field value a; s |A|, correspond- 
ing to the magnetic-ion spin splitting. 

The branch of stiff excitations f^ iff , Eqs. Q28JI and 
(|31 | , is primarily band- like in character and is centered 
around the much larger energy scale of the itinerant car- 
rier mean-field spin splitting |A|. 



D. Moderate Mn doping 

We now switch to the case of higher dimensions 
(N + ,N~ > 1), which corresponds to several Mn ions 
across the quantum well. As in the previous Sec. If 11 CI 
dedicated to dilute Mn doping, we consider Mn spins cou- 
pled ferromagnetically to the itinerant carriers (A < 0) 
at T = 0. Comments regarding the antifcrromagnetic 
coupling case will be introduced later in Sec. IIVI As 
an example we choose N + = N~ = 5 with a relative 
spin density x s = 0.05. The dispersions corresponding to 
even and odd modes are obtained numerically by solving 
det[D- n 1 / (k,n)]+ 8 = and det^fk.O)]^ = 0, re- 
spectively. This leads to a set of six even modes, five rela- 
tively soft (flm) and one hard (fi^ iff ), and five odd modes 

(f2^). The index I (1 < I < 5) orders the modes from 

the bottom (I = 1) to the top (I = 5) of the spectrum 
in each case. Our results for £1%. (solid lines) and fiZ. 
(dashed lines) are summarized in Fig. Eta)-(c) for three 
characteristic ratios e s /|A| = 0.9, 0.975, and 1.5, respec- 
tively. Panels (a) and (b) correspond to half-metallic car- 
riers while panel (c) depicts results for the case of partly 
polarized carriers. Related results for iU ti a are shown in 
Fig. El In all plots the shaded zones represents the Stoner 
continuum. For each case, the normalized excitation en- 
ergies fV|A| are plotted as a function of the normalized 
k-vector ^/e k /|A| 



respond to k/k-o <C 1, the regime well below the Debye 
cut-off in which the continuum approximation is most 
reliable. 



0.05 
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0.05 
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(b) 




A(k/k D ), where A = y/^J\A\. 
For typical DMS quantum well systems we estimate that 
A > 10. Hence, the results shown in Figs. and 0] cor- 



/£kW 



FIG. 3: Dispersion of low-energy spin excitations in a single- 
subband DMS quantum well with ferromagnetic coupling 
(A < 0) and dimensional cut-off N + = iV~ = 5 (see text). 
Panels (a) and (b) correspond to half-metallic (i.e. fully po- 
larized) carriers (e a /|A| = 0.9 and 0.975, respectively). Panel 
(c) depicts results for partly polarized carriers (e s /|A| = 1.5). 
The even modes are denoted by solid lines and the disper- 
sionless odd modes by dashed lines. The spatial profiles 
of the even and odd modes are illustrated in Fig.^Ja)-(e) and 
Fig. |7| respectively. The integer labels order the modes by 
increasing frequency. The shaded zones indicate the Stoner 
continuum. 

We start by discussing the properties of the relatively 
soft even modes Slt^ depicted in Fig. |3 by solid lines. 
Panel (a) corresponding to e s /|A| = 0.9 is representa- 
tive of results for half-metallic cases with e s /|A| <C 1. 
There we find a set of modes at k = that are dis- 
tributed in the small energy window < Q < 2x s \A\ and 
have k 2 dispersion at finite wavevectors, corresponding 
to a finite spin-stiffness. The upper limit in this spec- 
trum is twice the mean-field value predicted in the case 
of bulk systems^ This is due to the fact tha t the car- 
rier spin density is modulated by x{ z ) = y/ 2 / d sisi(tt z / d) 
across the quantum well (lower density close to the border 
z = 0, d and higher density close to the center z = d/2). 
Spin modes with large relative amplitude near the center 



FIG. 4: Dispersion of high-energy spin excitations in a single- 
subband DMS quantum well with ferromagnetic coupling 
(A < 0) and dimensional cut-off N + = 5 (see text) for half- 
metallic (i.e. fully polarized) carriers with e s /|A| = 0.9. Sim- 
ilar features appear for partly polarized carriers (e s /|A| > 1). 
The stiff branch fi+ ;ff (solid lines) corresponds to the mode 
w stifi (Fig-IMf))- The Stoner continuum is represented by the 
shaded zone. Ferromagnetic interactions with the local mo- 
ments peel a collective mode off the particle-hole continuum. 



of the quantum well see a carrier spin density which is 
effectively higher than in the uniform spin-density bulk 
case. 

The first branch £27L corresponds to a gapless 
Goldstone-mode u/^ (z) = u>o(z), Fig-E^a), similar to the 
one discussed above for the dilute Mn doping case, in Sec. 
If II CI For the I > 1 modes, the spatial structure (z) 
is not constant within the quantum well (see Figs. EJTb)- 
(e)) and the energies O^(k) therefore approach a finite 
value as |k| — > 0. The dependence of the excitation en- 
ergies n^(k = 0) on I is not obvious, since the gap is 
not simply related to the effective transverse momentum 
kf^ = j Q dz U3%,dw%./dz associated to each mode. In- 
stead, the local excitation density w~^{z) and its cor- 
relation with the carrier density, proportional to x 2 ( z )j 
is more relevant. (See Fig. [S] for a comparison). We 
illustrate the i-dependence of fi^(k = 0) in Fig. for 
N + = 20 (full circles). We observe a spin-mode accu- 
mulation close to the boundaries of the spectrum that 
is not evident for small N + . Furthermore, the inset in 
Fig. shows that the dispersion does not appear to be 
quadratic for small I. Since f2^(k = 0) = 0, a quadratic 
fit has a single free parameter. The dashed (dotted) line 
corresponds to a quadratic fit between the I = 1-mode 
and the I = 2- (I = 5) one. The curves differ by a (rela- 
tively large) factor of order 3.5. 

As e s /|A| increases (Fig. Etb)) and the Stoner contin- 
uum meets the different f2^ branches, the correspond- 
ing spin stiffness drops to zero. For e s /|A| > 1 (partly 
spin-polarized carriers, Fig.^c)) all branches are nearly 
dispersionless until they enter the particle-hole contin- 
uum; the softness of these excitations is certain to have 
an impact on the magnetic properties of these ferromag- 
nets. The spectral density shows that the relatively soft 
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FIG. 5: Spatial profiles across a single-subband DMS quan- 
tum well for even parity collective modes (solid lines). The 
well has width d (0 < z/d < 1), ferromagnetic coupling 
(A < 0), and the dimensional cut-off N + — 5 (see text). The 
carrier density X 2 ( z ) (dashed line) is shown for comparison. 
Note that the Goldstone mode, panel (a) , is constant in space 
corresponding to uniform spin rotation. The other modes 
tend to have higher energy when they have higher weight 
toward the middle of the quantum well, where exchange in- 
teractions are stronger. 
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FIG. 6: Dispersion of low-energy spin excitations in a single- 
subband DMS quantum well for in-plane momentum k = 0, 
ferromagnetic coupling (A < 0), and dimensional cut-off 
N + — N~ — 20, as a function of the mode index I (see text). 
The full circles depict the excitation energies Qf„ (k = 0) 

of even modes while the empty circles correspond to 

f2^(k = 0) for odd modes oCs. The inset corresponds to 
two extreme quadratic fittings (dashed and dotted curves, see 
text) to the dispersion of even modes (full circle) for small I. 
The curves differ by a prefactor of order 3.5. 



modes survive in the midst of the Stoner continuum, but 
the stiff mode does not as discussed below. 

Regarding Mn doping and kernel dimensionality, N + 
increases with higher magnetic-ion density and quantum 
well width. In this situation, new, relatively soft, modes 
arise from the top of the spectrum squeezing the rest to 
the bottom in order to satisfy < fn\ < 2x s |A|. This 
happens because increasing the dimension of the kernel 
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admits the presence of higher order Fourier components 
in our expansion and allows lower energy modes to aquire 
a larger relative amplitude at the borders of the quantum 
well, where the free-carrier density is reduced and the Mn 
spin splitting is consequentially weaker. The same con- 
sideration applies to the odd mode (see below) behavior 
as a function of N ~ . 

In addition to the relatively soft modes we find one 
even, stiff branch iU ti a (Fig. lying above the Stoner 
continuum. Its properties are similar to those discussed 
for the stiff mode in the dilute Mn doping case, Sec. IIII Cl 
The mode is restricted to relatively large wavelengths if 
compared with the case of soft modes. This is due to 
the proximity of the Stoner continuum and their strong 
interaction. The corresponding excitation profile w^ iff (,z) 
is shown in Fig. [SJf ) (solid line) . Its similarity to the 
carrier density profile x 2 ( z ) (dashed line) demonstrates 
that this mode is primarily associated with the itinerant- 
carrier subsystem dynamics. 

We continue discussing briefly the properties of the 
dispersionless odd modes Q7* depicted in Fig. [3] (dashed 
lines). As pointed out above, the flat dispersions have 
their origin in the absence of correlation effects related 
to spin reorientations. Odd mode fluctuations of the lo- 
cal moment give rise to effective fields that are averaged 
to zero by x 2 (z) and are consequently not correlated with 
carrier spin fluctuations. This makes the odd modes 
transparent to the Stoner excitations (they do not in- 
teract with the particle-hole continuum, unlike the even 
modes) and independent of the ratio e s /|A| (that is, of 
the carrier spin polarization). Moreover, all modes are 
lodged within the (low) energy window < ft < 2a; s |A| 
and present a finite gap whose magnitude depends on 
the particular excitation profile u^(z) (Fig. 0). This 
behavior is similar to that found for the relatively soft 
even modes except that there is no Goldstone mode in 
the odd spectrum; the higher the weight at the border 
of the quantum well, the lower the gap (see e.g. uj^Az), 

Fig.0e), which corresponds to in Fig.|3). As in the 
even case, the dependence of 0,7* on I is nontrivial. See 
Fig. HJ1 (empty circles) for an illustration of this depen- 
dence in the case N~ = 20. 

We further note that, as can be seen from Eq. (|20|l . in 
the limit of a large number of vertical modes (N + , N~ 3> 
1) it holds Q7-. w fii +1 x(k = 0), provided that x s <C 1. 
The even and odd mode pairs are then nearly degener- 
ate even and odd combinations of excitations at opposite 
edges of the quantum welli^i 

Our analysis holds for symmetric quantum wells only. 
We comment shortly on how an asymmetry of the 
quantum-well potential will affect or results. On the one 
hand, a gapless Goldstone mode as found for symmetric 
quantum wells when neglecting spin-orbit coupling (see 
also Sec. IIII E|l still exists. On the other hand, the spin- 
wave modes will have no definite parity anymore, and 
the above classification into even (dispersive) and odd 
(nondispersive) modes no longer holds. The particular 



excitation profiles (Figs. |3]and[?|) and energies (Fig. 0) 
will depend on the local-carrier density \ 2 { z )- 
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FIG. 7: Spatial profiles for spin excitation modes with odd 
parity (solid lines) across a single-subband DMS quantum well 
of width d (0 < z/d < 1) with ferromagnetic coupling (A < 0) 
and dimensional cut-off N~ = 5 (see text). The transverse 
carriers density X 2 ( z ) (dashed line) is shown for comparison. 



E. Effect of spin-orbit coupling 

The presence of spin-orbit coupling described by the 
Dresselhaus Hamiltonian H^, Eq. JSJ, introduces easy- 
axis magnetic anisotropy which is, of course, necessary 
for long-range magnetic order in a quantum well. When 
the anisotropy, which explicitly breaks rotational invari- 
ance for the magnetization orientation, is accounted for 
a finite energy-gap f2 so appears in the lowest lying collec- 
tive mode branch and several of the lower-lying branches 
in Figs. El and are shifted to higher energies. We calcu- 
late 51 so for the lowest even mode of constant excitation 
profile wt-., Fig.Efa). With this aim we follow the proce- 
dure of Sec. IIII 51 and calculate fi/jj for small spin-orbit 
coupling and k = 0. We find that 



a 



(32) 



SN Mn d max{|A|,2e F }' 
where ep = (ti 2 /2m*)kp is the Fermi energy of the 2D 



electron gas paramagnetic state and n 



iK This 



coincides with our previous resulli^ obtained by using 
perturbation theory. 

The above result and discussion holds for symmetric 
quantum wells. In asymmetric quantum wells, spin-orbit 
coupling is not only described by the Dresselhaus terms, 
but has also a RashbaS spin-orbit contribution. For elec- 
trons in the conduction band, the Rashba Hamiltonian is 
linear in the in-plane momentum k. This will affect the 
magnetic anisotropy and the energy gap in the spectrum 
of collective spin excitations. The interplay of these two 
types of spin-orbit coupling will complicate the disper- 
sion of the spin waves (Fig. |5J). In particular, they will 
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become anisotropic, analogous to the anisotropic trans- 
port properties discussed in Ref. For heavy holes 
in the valence band the leading Rashba term is cubic 
in momentum^ This complicates the evaluation of the 
spin-wave dispersions even more. 



IV. COMMENTS AND CONCLUSION 

We have developed a theory of collective spin excita- 
tions in DMS quantum wells by extending the approach 
that we used previously for bulk systems «2£ The theory 
goes beyond mean-field and RKKY approaches and ac- 
counts for both finite itinerant carrier spin-splitting and 
dynamical correlations. We applied this tool to the study 
of spin excitations in the ordered magnetic state at zero 
temperature. As in the bulk case, we have recognized two 
different energy scales on which the spin excitation spec- 
trum depends: one hard scale |A| principally related to 
the itinerant-carrier subsystem, and one soft scale ir s |A| 
for the magnetic- ion spin excitations, where x s |A| <C |A|. 
In addition, a continuum of Stoner excitations (corre- 
sponding to flipping a single spin in the itinerant-carrier 
band) also emerges from the theory. Although most rele- 
vant to DMS ferromagnet properties in circumstances for 
which the magnetic moments have a high degree of spin 
alignment, this theory of the elementary excitation spec- 
trum of the system sheds considerable light on the nature 
of the magnetic state and on the physics that controls the 
critical temperature of the system. 

The excitation spectrum of this magnetic system is 
quite unusual because of its ambiguous dimensionality. 
A slab of magnetic ions is coupled by a 2D electron sys- 
tem that is frozen into a single growth direction electronic 
subband and cannot distort its z-dependence to accom- 
modate magnetic fluctuations. We find that the excita- 
tion spectrum of this system has multiple 2D branches. 
The number of reasonably well-defined branches of ex- 
citations that have primarily local moment character is 
close to the width of the quantum well measured in units 
of the mean-separation between Mn ions, as expected by 
analogy with a reference systems in which the local mo- 
ments are placed on a lattice with the same volume per 
Mn and a finite number of layers. On the other hand, we 
find that there is only one 2D branch of collective exci- 
tations that have primarily electronic character. When 
spin-orbit interactions are neglected, the gapless Gold- 
stone mode branch has quartic rather than quadratic dis- 
persion, implying that the spin-stiffness vanishes, except 
when the carrier system is half-metallic. Unless, that is, 
e s < | A | and the mean- field carrier spins are consequently 
fully spin-polarized. This property arises somewhat ac- 
cidentally from the particular features of effective inter- 
actions mediated by carriers in 2D parabolic bands and 
has partly been noted for the RKKY limit (e s ^S> |A|) in 
previous work^ If we had included carrier-carrier inter- 
actions in our theory, the spin-stiffness would not vanish 
but, depending on the carrier density, might have a neg- 



ative sign implying that the ferromagnetically ordered 
state is unstable. When spin-orbit interactions are in- 
cluded, however, the collective excitation spectrum has a 
small gap and negative dispersion in the lowest-lying col- 
lective mode, while unusual, does not necessarily imply 
instability. Finally, we remark that the excitation spec- 
trum includes a large number of even and odd weakly dis- 
persive or non-dispersive branches in which fluctuations 
in the local moments are concentrated in particular parts 
of the quantum well and do not couple strongly to fluc- 
tuations in the band-electron spin-orientation. In these 
modes, the excitation energy is determined primarily by 
the local strength of the mean-field interaction between 
the Mn moments and the band electrons, which becomes 
small because of carrier quantum-size effects toward the 
edges of the quantum well. 

The results summarized above for the collective exci- 
tation spectrum suggest that thermodynamic properties, 
the temperature dependence of the magnetization for ex- 
ample, are likely to be quite unusual in DMS ferromag- 
nets, particularly when the carriers are in the conduction 
band where spin-orbit interactions are rather weak. The 
influence of thermal fluctuations on the magnetization 
will be enhanced not only by the reduced dimension 23 , 
but also by the small and possibly even negative spin stiff- 
ness mentioned above. (Fluctuations in the Goldstone 
mode branch have a relative importance that goes like 
X/N , where N is the effective number of Mn layers in the 
film that we have discussed previously.) In the case of va- 
lence band DMS quantum well fcrromagnets, strong mag- 
netic anisotropy 55 should lead to ferromagnetism that is 
essentially Ising in character. The magnetization should 
then be fairly constant over a wide interval of tempera- 
ture before dropping fairly rapidly to zero near the crit- 
ical temperature. In the conduction band case for which 
the model we have studied applies most directly, how- 
ever, the gap in the excitation spectrum will be quite 
small, much smaller than the mean-field ferromagnetic 
critical temperature, as we have discussed earlier^ The 
true critical temperature is likely to be determined in 
large measure by long length scale fluctuations and to be 
substantially smaller than the mean-field temperature. 
Since the stiffness of this system is very small, it will be 
significantly altered by spin-orbit interactions; this part 
of the physics is something that we have not addressed 
here. 

Furthermore, the model we have studied ignores dis- 
order, which is likely to play an important role in adju- 
dicating the way in which these subtle competitions are 
resolved. We believe that careful study of the magne- 
tization and other properties of electron quantum well 
DMS ferromagnets at low temperatures will reveal a lot 
of subtle, unusual, and interesting physics. 

Finally we emphasize that the numerical results pre- 
sented here are for the case of ferromagnetic interactions 
between the carriers and the local moments (A < 0), 
which is expected to apply for n-doped semiconductors. 
Since we have taken the local-moment spins to point 
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down in their ground state, this means that the up spins 
are the minority spins and the down spins are the ma- 
jority spins in the ferromagnetic case while their roles 
are interchanged in the antiferromagnetic case. It fol- 
lows that the only change in Eq.^Jwhen A changes sign 
is that iv n changes sign. This change has a number of 
consequences that are fairly subtle when x s is small, but 
can in principle be more consequential. The most impor- 
tant differences are that the collective mode with dom- 
inant carrier character which appears above the Stoner 
particle-hole continuum in the case of ferromagnetic in- 
teractions, lies below the Stoner continuum in the case of 
antiferromagnetic interactions. In addition 1/ (l+x s ) fac- 
tors which appear in expressions for the collective mode 
energies, factors that express the band electron contribu- 
tion to the total spin density of the system, are replaced 
by 1/(1 — x s ) factors in the antiferromagnetic case. These 
factors are not present in a RKKY description of the car- 
rier mediated interactions. 
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APPENDIX A: MEAN-FIELD ITINERANT 
CARRIER SPIN POLARIZATION 

The 2D mean-field itinerant carrier spin density for 
spin a and subband spin-splitting A = J cx NmiiS in 



the presence of weak Dresselhaus spin-orbit coupling 
at B = is given by n 17 = j d 2 k / '(2tt) 2 / (e£) , where 
€ = £k - (<t/2)A + cre k ° - //, e k = h 2 k 2 /(2m*), e k ° = 
— 7 2 (fc 2 ) 2 fc 2 /A, p,' is the effective chemical potential of 
the 2D carrier gas, and /(e) is the Fermi distribution. 
For zero temperature we find 



where we have used /(e) t=o = $( — e )- The difference 
between up and down contributions determines the net 
carrier spin density defined as p = (r?J — n^)/(n) -\-<nr). 
Denoting the Fermi energy of the majority-spin carrier 
band by e s = p! + |A|/2, we see from Eq. (|A1|I that the 
carrier system is fully spin-polarized or half-metallic (i.e. 
\p\ = 1) when e s < |A|. On the other hand, partly spin- 
polarized carriers (i.e. \p\ < 1) correspond to e s > |A|. 
APPENDIX B: CORRECTION DUE TO 
CORRELATION EFFECTS 



Correlation effects due to the response of quantum 
well carriers to Mn spin reorientations are taken into ac- 
count in the kernel Z? -1 (k, z, z' , ft) by the second term of 
Eq. I|14f) . This is manifested by the momentum and en- 
ergy dependence of the integral /(k, ft), Eq. I|23|) . In the 
absence of magnetic field and spin-orbit coupling we find 
for zero temperature (see Appendix lAl for definitions) 



d 2 9 [/(4) - /(< +k )] T =o 



(27T) 



n + 4 - <+k 



4irh 2 e k 



6(p'-A/2) (ft + A-e k ) l-Wl 



4(/i' - A/2)e k 
(ft + A-e k ) 2 



(V + A/2) (ft + A + e k ) 1-4/1 



4(/i' + A/2)e k 



(fi + A + e k ) 



,1 



n l (p! + A/2 + n) - n T (p! - A/2 - ft) 



ft + A (ft + A) 3 

-(nT+n i )/e k + 0(l/e 2 ), 

I 



e k + 0(e 2 ) 



(Bl) 

(B2) 
(B3) 



where the integral has to be considered as a Cauchy prin- 
cipal value. Eqs. (|B2|I and l|B3l) correspond to first or- 
der expansions in e k and l/e k , respectively. The prefac- 



tors that express the ^-dependence in the second term of 
Eq. (|14l) contain information on the quantum well geom- 
etry and lead to M 2 in Eq. J5UJ). 
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